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ABSTRACT 


A  proRrjua  has  been  written  for  the  HP-97  (HP-67)  minicomputer  to  solve 

the  blast  wave  from  a  nuclear  detonation.  The  program  first  determines 

standard-altitude  paraoMtere  for  the  altitude  of  interest.  Using  these 

paraoietera  as  input,  the  program  then  calculates  pressures,  temperatures, 

densities,  velocities,  and  Mach  Numbers  in  front  and  behind  the  shock. 

C<  or  s 


Three  subroutines  allows  the  user  to  input  (for  any  altitud^i|^2,000  ft)  a 
specific  overpressure  across  the  shock,  a  specific  gust  velocity  behind  the 
blast  wave,  or  a  specific  dynamic  pressure  behind  the  blast  wave.  Corre¬ 
sponding  parameters  for  each  input  are  calculated. 


NUCLEAR  BLAST  PROGRAM 


For  Mini-Calculators 

Introduction 

In  the  process  of  perforwino;  calculations  of  the  values  of  the  pres¬ 
sures,  temperatures,  velocities,  densities,  etc.  behind  weak  shock  waves  at 
some  diatfinco  f  i  c  m  a  nuclear  detonation,  it  became  obvious  that  hand  calcu¬ 
lation  was  tedious  and  interpolation  from  {graphical  results  subject  to 
e-ror,  Therefot-ft  a  simple  prop, ram  was  developed  for  use  on  the  HP-97  (or 
HP-i7)  mini  calculator  to  perform  these  operations.  Ihis  v>ropram  is 
describftci  in  the  lollowinp  narrative. 

Although  this  simple  program  is  not  intended  to  replace  the  more 
sophisticated  contpoter  coder;  o-rHll.ible  for  such  computation,  it  is  a 
valuable  tool  for  first  cvit  analyses.  It  also  h-tis  the  distinct  advantage  of 
being  highly  portable,  and  serves  as  bandy  reference  for  use  when  the 
detailed  data  are  not  immediately  available. 

Discussion 

At  distances  from  any  nuclear  detonation  of  interest  to  the  analyst 
working  with  airborne  systems,  the  blast  wave  is  relatively  weak  and 
amenable  to  solution  via  standard  fluid  mechanical  equations.  The  first 
step  in  such  s  solution  is  tVie  application  of  standard  fluid  dynamical 
equations  used  to  solve  the  stationary  shock.  See  Figure  1.  The  governing 
equations  relating  conditions  in  front  of  and  behind  Che  shock  can  be 
obtained  from  any  text  on  supersonic  flow.  I  used  Shapiro*. 


*  A.  H.  Shapiro,  The  Dynamics  and  Theraodynaaica  of  Compresaible  Fluid  Flow. 
Honald  Presm,  N.Y.  1953 


iA,>i-V>n\. 


Figure  1. Flow  parameters  Cor  a 
stationary  shock  wave 


After  solving  the  stationary  shock  problem,  1  then  transformed  the 
problem  into  a  moving  shock,  or  blast  wave  problem  by  allowing  the  shock  to 
move  with  a  speed  such  that  its  Mach  number  was  equal  to  •  This  system 
uses  primed  variables  and  is  depicted  in  Figure  2. 
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I  Figure  2.  Flow  parameters  for  moving  shock  in 
terms  of  stationary  shock  solution. 


As  depicted  in  Figure  2,  the  stationary  shock  solution  provides  all  the 
data  needed  to  solve  the  blast  wave  problem.  The  parameters  of  most 
interest  to  the  nuclear  survivability/vulneiability  analyst  are  V2'  (the 
gust  velocity),  the  overpressure,  (P2*  ”  '  “  ^2  ” 

dynamic  pressure,  q2'. 
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Tlie  ptogrart!  incorporating  the  solution  ot  ttie  blast  wave  is  listed  in 
the  Appendix  H  and  detailed  operating  instructions  are  contained  in  Appendix 
A,  The  ''gutr.”  of  the  program  is  rubroutioe  B,  which  solves  the  stationary 
shock  then  utilizes  that  soliition  to  solve  the  blast  wave  (moving  shock)  for 
any  selected  altitude  and  overpressure.  Subroutines  C  and  D  accept  gust 
velocity  and  dynamic  pressure  inputs  respectively  and  iterate  back  through 
the  V'Qsio  part  o'  ?uj:i' tontine  B  until  the  calculated  quantity  is  within  the 
acceptable  tolerance  of  the  input. 

Subroutine  A  is  the  cornoretone  of  the  program.  Tt  confains  analytical 
expressions  describing  the  variation  of  the  standard  atmosphere  with 
altitude  as  presented  in  Dommasch.*  For  the  specified  altitude  these 
equal 'ons  are  solved,  resulting  in  ratios  of  the  values  of  the  parameter  of 
interest  at  aUitude  to  the  level  standard  value.  Multiplication  of 
these  ratios  by  the  pertinent  sea  level  standard  results  in  the  desired 
value  at  altitude.  These  values  are  stored  and  also  used  in  subsequent 
calculations.  Note  that  because  of  program  step  limitations,  the  program  was 
separated  into  two  parts.  Patt  1  it'  po'i.inent  for  the  troposphere  (b<^ 
36,000  feet)  whi  le  part  2  i  i-  aocMrat-:  in  th«»  stratosphere  (36,000<h 
f}2,000)  wliere  the  temperature  is  constant  with  altitude. 


CONCLUSION 

The  program  presented  herein  should  prove  to  be  a  valuable  tool  to  analysts 
dealing  with  weak  blast  waves  of  any  type,  but  should  prove  of  particular 
interest  to  nuclear  survivability/vulnerability  analysts  in  government  and 
Industry. 


*D.  0.  Dommasch,  S.S.  Sherby,  and  T.F.  Connolly,  Airplane  Aerodynami< 
Pitman  Publishing  Corp.  N.Y.  1961. 
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I.  Subroutine  A,  Calculation  of  Aobient  Ataoepheric  Paraaetera. 

Upon  loading  the  program,  the  first  step  must  be  the  utilization  of 
Subroutine  A  to  calculate  the  ambient  static  pressure  and  temperature  at  the 
altitude  of  interest.  The  other  subroutines  use  these  parameters,  and  will 
not  function  unless  the  calculations  are  made.  Note:  Load  Part  1  if 
altitude  of  interest  is  36,000  feet  or  lees,  load  Part  2  if  the  altitude 
of  interest  is  between  36,000  and  82,000  feet. 

PRESS  A  -  This  step  clears  all  registers  and  sets  up  the  program. 

ENTER  ALTITUDE  -  Altitude  is  in  feet. 

PRESS  R/S  -  Calculation  requires  a  few  seconds.  When  the  altitude  is 
displayed,  calculations  are  complete.  Secondary  registers  contain  results 
of  the  calculations. 
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To  recall  ihese  pgratn . t ers ,  press  t,  then  P^S.  Use  RCL  for  register  of 
Interest.  Pi .  to  uho  of  other  aubroutines  either  restore  the  original 
configuration  (I,  P^S)  or  restart  program. 

II.  Subroutine  B  -  Calculation  of  Gust  Velocity  (given  an  Overpressure). 
This  suhroutine  accepts  an  overpressure  (in  psi)  and  solves  the  blast 

wave.  The  gusi  vlnciry  is  displayed  at  the  completion  of  the  calculations. 
(Recall  that  Subroutine  A  must  be  exercised  at  least  once  prior  to  use  of 
this  subroutin**.  The  results  of  the  subroutine  are  pertinenf  only  for  the 
last  altitude  entered  in  Subroutine  A.) 

PRESS  R  Initializes  subroutine 

ENTER  OVERPREilSllRK  In  psi.  Recall  overpressure  “  ^2  ”  ^1* 

PRESS  R/S  Arler  ca J culaf.ions  are  complete,  gust  velocity  is 

d i  spl  .syed  . 

At  the  point  the  user  has  several  options.  If  the  gust  velocity  corre¬ 
sponding  to  the  selected  overpressure  for  a  different  altitude  is  desired, 
press  A,  enter  the  nt;w  altitude,  prea.s  R/S,  then  exercise  subroutine  B 
again.  If  Lous  over HU'esf-:*!  e/gust  v.-ilocity  pairs  are  desired  for  the 
selected  altitude,  it  is  not  necessary  to  recalculate  the  atmospheric  param¬ 
eters,  press  B,  then  enter  the  overpressure.  The  corresponding  gust  is 
displayed . 

III.  Subroutine  C.  Calculation  of  Overpressure  (given  a  gust  velocity). 
This  subroutine  accepts  a  gust  velocity  (in  feet/second)  and  calculates 

the  corresponding  overpressure  via  iteration  through  the  basic  solution 
incorporated  in  subroutine  B.  This  subroutine  may  be  exercised  immediately 
after  the  altitude  computations,  or  the  completion  of  any  of  the  other 


subroutines . 


PRESS  C 


Initializes  subroutine 


ENTER  GUST  VELOCITY  feet/second 

PRESS  R/S  Several  iterations  may  be  necessary  for 

convergence.  After  the  calculations  are  complete, 
the  overpressure  is  displayed. 

If  other  gust  velocity/overpressure  pairs  are  desired  for  the  same  altitude, 
simply  repeat  the  above  steps.  If  a  new  altitude  is  desired,  subroutine  A 
must  be  reaccomplished  prior  to  initialization  of  subroutine  C. 

IV.  Subroutine  D.  Calculation  of  the  Overpreaaurc  and  Guat  Velocity 
(given  a  dynamic  presaure). 

This  subroutine  accepts  a  dynamic  pressure  in  psi  and  calculates  the 
corresponding  overpressure  and  gust  velocity  via  interation.  This 
subroutine  may  be  exercised  iramedirtely  .after  the  completion  of  subroutine 
A,  or  upon  completion  of  any  of  the  other  suhroutines. 

PRESS  I)  Initializes  subroutine, 

ENTER  DYNAMIC  PRESSURE  in  psi.  NOTE.  A  maximum  of  3  places  after  the 

decimal  will  be  tolerated  -  otherwise  program 

cycles  endlessly. 

PRESS  R/S  Several  iterations  may  be  required  for 

convergence  after  calculations  are  complete, 
overpressure  is  displayed, 

PRESS  R/S  Gust  velocity  is  displayed. 

V.  Miscellaneous  Notes 

1.  If  the  program  cycles  endlessly,  or  if  an  error  is  indicated  in  the 
display,  press  A  twice  to  clear. 
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Appendix  B 
PROGRAM  LISTING 


I  Part 

001 

1 

*LBLA 

21  11 

002 

RCL9 

36  09 

003 

CLRG 

16-53 

004 

P  S 

16-51 

005 

CLRG 

16-53 

006 

R/S 

51 

007 

ST09 

35  09 

008 

6 

06 

009 

• 

-62 

010 

8 

08 

011 

7 

07 

012 

5 

05 

013 

EEX 

-23 

014 

6 

06 

015 

CHS 

-22 

016 

X 

-35 

017 

CHS 

-22 

018 

1 

01 

019 

+ 

-55 

020 

STO0 

35  00 

021 

4 

04 

022 

• 

-62 

023 

2 

02 

024 

5 

05 

025 

6 

06 

026 

1 

01 

027 

YX 

31 

028 

ST02 

35  02 

029 

ROL0 

36  00 

030 

5 

05 

031 

• 

-62 

022 

2 

02 

033 

5 

05 

034 

6 

06 

035 

1 

01 

036 

YX 

31 

037 

STOl 

35  01 

038 

RCL0 

36  00 

039 

5 

05 

040 

1 

01 

041 

8 

08 

042 

X 

-35 

043 

ST03 

?5  03 

044 

RGLl 

36  01 

045 

1 

01 

046 

4 

04 

047 

• 

-62 

048 

7 

07 

049 

X 

-35 

050 

ST04 

35  04 

051 

RCL2 

36  02 

052 

* 

-62 

053 

0 

00 

054 

0 

00 

055 

2 

02 

056 

3 

03 

057 

7 

07 

058 

8 

08 

059 

X 

-35 

060 

ST05 

35  05 

061 

RCL3 

36  03 

062 

X 

54 

063 

4 

04 

064 

9 

09 

065 

X 

-35 

066 

ST06 

35  06 

067 

RCL3 

36  03 

068 

STOC 

35  13 

069 

RCL4 

36  04 

070 

STOA 

35  11 

071 

RCL9 

36  09 

072 

P  S 

16-51 

073 

ST09 

35  09 

074 

*LBLB 

21  12 

075 

R/S 

51 

076 

STOB 

35  12 

077 

*LBL2 

21  02 

078 

RCLA 

36  11 

- - 

- - 

--  ■ 

079 

-•55 

108 

-62 

137 

X 

-35 

080 

STOD 

35  14 

109 

2 

02 

138 

RCLE 

36  15 

081 

KCLA 

36  11 

110 

X 

-35 

139 

X 

54 

082 

•r 

-24 

111 

1 

01 

140 

- 

-45 

083 

ST02 

35  02 

112 

4- 

-55 

141 

ST06 

35  06 

084 

2 

02 

113 

STOl 

35  01 

142 

RCLC 

36  13 

085 

t 

-.62 

114. 

RCLD 

36  14 

143 

RCLl 

36  01 

086 

4 

04 

1  15 

• 

-62 

’44 

X 

-35 

087 

X 

-35 

116 

2 

02 

145 

X 

54 

088 

• 

•  62 

117 

X 

-35 

146 

4 

04 

089 

4 

04 

118 

1 

01 

147 

9 

09 

090 

+ 

-55 

1  19 

4 

-55 

148 

X 

-35 

091 

2 

02 

120 

RCLl 

36  01 

149 

RCL6 

36  06 

092 

t 

-62 

121 

-24 

150 

X 

-35 

093 

8 

08 

122 

STOl 

35  01 

151 

ST07 

35  07 

094 

A 

-24 

123 

RCL2 

36  02 

152 

RCLA 

36  11 

095 

STOD 

33  14 

124 

RCLl 

36  01 

153 

RCLB 

36  12 

096 

5 

'  05 

125 

-24 

154 

+ 

-55 

097 

■f 

-55 

126 

ST03 

35  03 

155 

ROL6 

36  06 

098 

.SI'OE 

35  15 

127 

RCLl 

36  01 

156 

X2 

53 

099 

7 

07 

128 

X 

54 

157 

X 

-35 

100 

RCLD 

36  14 

129 

ST04 

35  04 

158 

• 

-62 

101 

X 

-35 

130 

RCL3 

36  03 

159 

7 

07 

102 

1 

01 

131 

1/X 

52 

160 

X 

-35 

103 

- 

-45 

132 

ST05 

35  05 

161 

ST08 

35  08 

104 

Rf 

36  15 

133 

RCL4 

36  04 

162 

RCL0 

36  00 

105 

A 

-24 

134 

1/X 

52 

163 

X-0? 

16-42 

106 

1/X 

52 

135 

RCLD 

36  14 

164 

GSBj  23 

16  15 

107 

STOE 

35  15 

13& 

X 

54 

165 

RCLl 

36  46 

12 


166 

X-0? 

16  42 

196 

RCLB 

36  12 

167 

GSB3 

23  03 

197 

*LBLD 

21  14 

163 

RCL7 

36  07 

198 

R/S 

51 

169 

*LRLC 

21  13 

199 

STOI 

35  46 

170 

R/S 

51 

200 

RCLB 

36  12 

171 

STO0 

35  00 

201 

X-0? 

16-43 

172 

RCL7 

36  07 

202 

2 

02 

173 

X-07 

16-42 

203 

STOB 

35  12 

174 

GSBe  : 

23 

16  15 

204 

GSB2 

23  02 

175 

2 

02 

205 

*LBL3 

21  03 

176 

STOB 

35  12 

206 

RCLI 

36  46 

177 

1 

01 

207 

RCLS 

36  08 

178 

0 

00 

208 

DSP3 

-63  03 

179 

2 

02 

209 

RND 

16  24 

180 

ST07 

35  07 

210 

X-Y? 

16-33 

181 

*LBLe 

21 

16  15 

211 

GSB5 

23  05 

182 

RCL0 

36  00 

212 

r 

-24 

183 

INT 

16  34 

213 

X 

54 

184 

RCL7 

36  07 

214 

RCLB 

36  12 

135 

INT 

16  34 

215 

X 

-35 

186 

X-Y? 

16-33 

216 

STOB 

35  12 

187 

GSBa 

23 

i  16  11 

217 

GSB2 

23  02 

188 

t 

-24 

218 

*LBL5 

21  05 

189 

RCLR 

36  12 

219 

0 

00 

190 

X 

-35 

220 

STOI 

35  46 

191 

STOR 

35  12 

221 

RCLB 

36  12 

192 

GSB2 

23  02 

222 

R/S 

51 

193 

*LBLa 

21  16  11 

223 

RCL7 

36  07 

194 

0 

00 

224 

R/S 

51 

195  STO0  35  00 
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